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ABSTRACT 

The evolution of X-ray emission from young massive star clusters is modeled, taking into 
account the emission from the stars as well as from the cluster wind. It is shown that the 
level and character of the soft (0.2-10 keV) X-ray emission change drastically with cluster 
age and are tightly linked with stellar evolution. Using the modern X-ray observations of 
massive stars we show that the correlation between bolometric and X-ray luminosity known 
for single O stars also holds for O+O and O+Wolf-Rayet (WR) binaries. The diffuse emission 
originates from the cluster wind heated by the kinetic energy of stellar winds and supernova 
explosions. To model the evolution of the cluster wind, the mass and energy yields from a 
population synthesis are used as input to a hydrodynamic model. It is shown that in a very 
young clusters the emission from the cluster wind is low. When the cluster evolves, WR stars 
are formed. Their strong stellar winds power an increasing X-ray emission of the cluster wind. 
Subsequent supernova explosions pump the level of diffuse emission even higher. Clusters 
at this evolutionary stage may have no X-ray bright stellar point sources, but a relatively 
high level of diffuse emission. A supernova remnant may become a dominant X-ray source, 
but only for a short time interval of a few thousand years. We retrieve and analyse Chandra 
and XMM-Newton observations of six massive star clusters located in the Large Magellanic 
Cloud. Our model reproduces the observed diffuse and point-source emission from these LMC 
clusters, as well as from the Galactic clusters Arches, Quintuplet and NGC 3603. 

Key words: stars: winds, outflows - stars: Wolf-Rayet - open clusters and associations: 
general - galaxies: LMC - X-rays: stars. 



1 INTRODUCTION 

Massive star clusters, and super star clusters (SSC) in particular, 
are among the most extreme cases of star formation regions in the 
local Universe. These clusters have a central stellar density up to 
~ 1O 5 M pc~ 3 and a total mass ranging from 10 3 to 10 6 M Q . 
Clusters of such mass and density can contain hundreds of thou- 
sands of OB stars, providing a unique laboratory to test the theory 
of evolution of massive stars. Moreover, these clusters can have a 
dramatic effect on their surrounding interstellar, and in some cases, 
intergalactic medium. Among Galactic examples of SSCs are the 
Arches and the Quintuplet cluster located in the vicinity of the 
Galactic center, and NGC 3603. A few dozens of massive star clus- 
ters are known in the Magellanic Clouds. 

Recent advantages in X-ray imaging and spectroscopy, result- 
ing in the detection of X-ray emission from some of the SSCs, 
prompted theoretical work to explain this emission. It is understood 
that regions of active star formation have high supernova (SN) 
rates. If the supernova energy input is th ermalised, a stro ng wind is 
driven out of the active region. IChevalier&ClegaNl985l) presented 
a self-similar analytical solution of the hydrodynamic equations for 
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the wind driven from a region of uniform mass and energy distribu- 
tion. This solution is scalable with the energy and mass input rates, 
and with the radius of the region of mass and energy production. 

Canto, Raga & Rodriguez (2000) considered a "cluster wind" 
resulting from the multiple interactions of stellar winds produced 
by the massive stars of a dense cluste r. Their cluster wind mode l 
is essentially the same as presented bv lChevalier & Clegd ll985l) . 
but accounts for mass loading due to stellar winds instead of super- 
novae. An analytical solution from of the hydrodynamic equations 
was obtained. A numerical simulation of a cluster of 30 stars was 
also presented. The authors applied their result to the Arches cluster 
and predict ed that the X-ray emission f rom its cluster wind can be 
detectable. lLaw & Yusef-Zader] 120041) have reported the discov- 
ery of diffuse X-ray emission fro m the Arches on a level similar as 
calculated bv lCanto et al.li200ol) . confirming their prediction. 

The evolution of the X-ray emission from a cluster of sin- 
gle young stars was investigated by Cervino, Mas-Hesse & Kunth 
(2002). They proposed an evolutionary synthesis model, where 
both supernova remnants and hot diffuse gas contribute to the X- 
ray emission. However, the hydrodynamics of the hot diffuse gas 
was not considered, and it was assumed that the X-ray luminosity 
of this gas is just some arbitrary fraction of the kinetic energy rate 
released in the cluster. 

A theoretical model which includes the effect of mass load- 
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Figure 1. XMM-Newton and Chandra images of star clusters in the LMC. The celestial positions of the clusters are indicated by arrows. R 136 is clearly 
visible in the top left frame, while the other five clusters are hard ly seen. Note that R 136, the brightest cluster, is by far the youngest cluster in the sample. The 
cluster coordinates and ages are from lMackev & Gilmord 120031) 



Table 1. X-ray observations of LMC massive star clusters 



Cluster 


Center (J2000.0) 


Detector 


Exp. time 


Count Rate 


name 


a 


<5 


name 


[ksec] 


[els" 1 ] 


R136 


05:38:42.5 


-69:06:03 


ACIS-I 


21 


8.31 x 10~ 2 


NGC2100 


05:42:08.6 


-69:12:44 


MOS 


6.5 


6' off-axis, non-det 


NGC 1818 


05:04:13.8 


-66:26:02 


ACIS-S 


25 


4.36 X 10~ 3 


NGC 1850 


05:08:41.2 


-68:45:31 


ACIS-S 


30 


on-axis, non-det 


NGC 1860 


05:10:38.9 


-68:45:12 


MOS1 


4.4 


off-axis, 5.9 X 10~ 3 


NGC 1831 


05:06:17.4 


-64:55:11 


ACIS-S 


25 


on-axis, non-det 



The coordinates of the cluster centers are from lMackev & Gilmord 120031) 



ing bas ed on the work of IChevalier & Clegd 1 19851) was consid- 
ered bv lStevens & Hartwelll 120031) . The authors compared the re- 
cent X-ray observations of local massive star clusters (Rosette, 
NGC346, NGC3603, R 136 and the Arches) with the predictions 
from the cluster wind theory. It was pointed out that from the ob- 
servational side the problems remain as to being sure that the dif- 
fuse emission is genuinely diffu se and associated with the cluster 
wind. lStevens & Hartwel}l2003l) showed that the diffuse X-ray lu- 
minosity of a cluster wind is correlated with the cluster wind ki- 
netic energy divided by the cluster radius, but the X-ray tempera- 
ture inferred from the available data is not well correlated with the 
predicted one. The authors concluded that from the available data it 
remains very unclear as to what is going on in the clusters regarding 
the X-ray emission. 

The research we are presenting here is stimulated by the X-ray 
images of the Galactic Center obtained with the XMM-Newton and 
Chandra observatories. Two rich massive star clusters, the Arches 



and the Quintuplet, have comparable masses, are located at approx- 
imately the same distance, suffer similar interstellar absorption, but 
their X-ray images, albeit being in the same field of view, are strik- 
ingly different. The Arches exhibits bright point sources and dif- 
fuse emission, while the Quintuplet is barely dete ctable. Only care- 
ful examination of the Quintuplet by [Law & Yusef-Zadebl 12004 
revealed the presence of weak X-rays from this cluster. There are 
hardly any point sources seen in the Quintuplet, and the ratio of dif- 
fuse X-ray to the IR flux is much higher in the Quintuplet than in 
the Arches. The X-rays from the Arches were modeled in a number 
of papers (see above), but the scarcity of X-rays from the somewhat 
older and less compact Quintuplet is not explained. 

Another sample of star clusters with similar structural param- 
eters, b ut different age is known in the Large Magellanic Cloud 
(LMC) iMackev & Gilmorel2003h . Conveniently, we could use the 
archival X-ray observations of six massive LMC clusters. Quite 
similar to the Galactic Center clusters, the youngest and most com- 
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pact cluster R 136 is bright in X-rays, while older clusters appear 
as non-detections, at least from the first glance. 

In the present paper we attempt to explain these observations 
as a consequence of cluster evolution. For ou r study we chose an 
approach similar to lStevens & HartwellH2003T) . but make use of the 
self-similar properties of the cluster wind solution. The effects of 
stellar evolution on the mass and energy deposition in a cluster are 
included via an evolutionary synthesis model. The important nov- 
elty is that we consider stellar X-rays as well as diffuse emission. It 
is often impossible to disentangle the diffuse and stellar emission in 
the dense cluster core. Therefore we estimate the luminosity of the 
stellar population and add it to the predicted luminosity of the clus- 
ter wind. The total, then, can be compared with observations. The 
stellar X-ray emission changes significantly when a star evolves, 
and we account for these changes as well. Our analysis is restricted 
to the inner parts of the clusters where the bulk of mass and en- 
ergy is deposited. This allows us to concentrate on the properties of 
cluster winds and stellar population, and to avoid considering the 
hot bubbles that are expected around regions with fast winds. 

The paper is organised as follows. The observations of mas- 
sive stellar clusters used in this work are described in Section 2. 
The approach to the modeling is specified in Section 3. The evo- 
lution of X-ray emission of cluster wind is calculated in Section 4. 
Section 5 is devoted to the evolution of the stellar X-ray luminosity. 
In this section we pay special attention to the correlations between 
X-ray and bolometric luminosities. In Section 6 we discuss the X- 
ray emission of supernova remnants, and conclusions are drawn in 
Section 7. 



2 OBSERVATIONS OF YOUNG MASSIVE CLUSTERS 

Recently a sample of six massive young star clusters located in 
the LMC was observed with the XMM-Newton and Chandra X-ray 
observatories. The observations are in the public archives (see Ta- 
blcQ. To our knowledge the results of these observations, except 
of R 136, have not been published. Therefore we retrieved the data 
on five LMC clusters from archives and analysed them. All clus- 
ters have similar masses, but differ significant ly in age. The struc- 
tural p arameters of the clusters are taken from iMackev & G ilmore 
120031) and listed in Tabled FigureQshows the X -ray images of the 
areas of the sky, where the clusters are located. Although being at 
the same distance and observed through similar absorbing columns, 
the clusters X-ray images are strikingly different. 

The youngest and most compact cluster, R 136, is at least one 
order of magnitude brighter than the more evolved clusters and 
displays wealth of bright point sources, but little diffuse emission. 
Among more evolutionary advanced clusters, three are not detected 
and the detected ones are apparently lacking point sources. 

In order to estimate observed X-ray fluxes from these LMC 
clusters, a standard analysis of pipeline processed products has 
been performed. Retrieved count rates were transformed to fluxes 
assuming that radiation comes from thermal bremsstrahlung with 
fcTx = 3keV. This assumption is based on the spectral fits of 
the diffuse emissi on from the Arches and th e Quintuplet clus- 
ters as discussed in Law & Yusef-Zadehli2004T) .The neutral hydro- 
gen column density for all LMC clusters was adopted as Nu = 
5 x 10 21 cm -2 , as derived from spectral fits of R 136 stars by Porte- 
gies Zwart, Pooley, & Lewin (2002). 

Below we briefly describe how the diffuse X-ray luminosity 
was estimated for each cluster under consideration. 

R 136 was observed by ACIS-I on board of the Chandra X- 



ray observatory for 21 ksec exposure time. However, even with the 
superb angular resolution of the ACIS camera (the pixel size is 
0.49"), it is hardly possible to separate stellar and diffuse emission 
within the cluster cor e, which has a radius 1.32" (corresponding to 
0.32 pc). However, as iMackev & Gilmord EobH) pointed out from 
the analysis of the cluster brightness profile, the inner region of the 
cluster is extended to 2.43 pc. We estimate the diffuse X-ray lumi- 
nosity of R 136 from this inner region. 

Using the wavedetect tool for the binned by factor of 6 
event file we obtained the total count rate in this region and 
translated it into the background subtracted flux, Fx = 1.6 x 
10 -12 erg s _1 cm ~ 2 . Background correction is important since 
R 136 is embedded in the extended region of diffuse X-ray emission 
30 Dor C. As the next step, we subtracted the lum inosities of the 
point sources CX1..9 IPortegies Zwart et alj|2002l) . and attributed 
the remaining luminosity to the diffuse emission. 

NGC2100 was observed 6' off-axis by XMM-Newton. Al- 
though clearly seen on an OM image, this rich young cluster is not 
detected in the X-rays (see Fig.0. An upper limit to the count rate 
from an area corresponding to the cluster core was obtained from 
a sensitivity map, and is 7 x 10~ 3 cts" 1 for the pn camera. Thus 
the upper limit on the flux from NGC2100 in the XMM-Newton 
passband is Fx < 3.5 x 10 -14 erg s _1 cm ~ 2 . 

The Chandra X-ray observatory has performed 25 ksec point- 
ing observations of NGC 1818 and NGC 1831, and a 30 ksec point- 
ing observation of NGC 1850. NGC 1818 is a faint source of ap- 
parently diffuse emission detected with the standard Chandra data 
analysis wavedetect tool. The estimated number of source counts is 
109. They are distributed nearly uniformly over a 15" x 15" area 
of the sky, with the center of diffuse emission coinciding with the 
cluster center as listed in the Table|2| Additionally, the psfratio de- 
termined by wavedetect is > 1, indicating that the source indeed 
may be extended. The low signal-to-noise spectrum of NGC 1818 
was extracted. 

In order to gain information on the absorption column density 
we extracted and fitted the spectrum of the brightest background 
object present in the ACIS-S field of view. The SeyfertI galaxy 
Cal F, only 10' aw ay from NGC 1818, is located behind the LMC 
iGeha et all2003l) . We tentatively fitted the spectrum of Cal F with 
a two-component model (power law and black body). This spectral 
fit allows to determine the absorption column density in the direc- 
tion to NGC 1818, N H « 1.4 x 10 21 citT 2 . One can make the 
plausible assumption that the absorption is mainly interstellar and 
is not intrinsic to the SeyfertI AGN. 

With the absorption column density constrained, the spectrum 
of the X-ray emission from NGC 1818 can be fitted either with a 
power law with index ~ 2.5, or with a multi-temperature thermal 
plasma with temperatures ranging from 0.3 to 5 keV. The raw and 
adaptively smoothed images of the cluster are shown in Fig.|2| The 
X-ray flux from NGC 1818 is Fx = 7 x 10~ 14 erg s" 1 cm ~ 2 . 

A similar observation of NGC 1831 and NGC 1850 failed to 
detect X-rays from these cluster. In order to estimate an upper limit 
on its X-ray flux we run the wavedetect tool with significantly lower 
detection threshold, therefore allowing for fake detections. Eventu- 
ally the area which coincides with the location of the clusters was 
"detected", providing an estimate of the count rate. The upper limit 
on the flux from NGC 1831 is F x < 4 x 10~ 14 erg s _1 cm ~ 2 , and 
from NGC 1850 Fx < 10" 14 erg s _1 cm " 2 . 

NGC 1860 was observed 9' off-axis by XMM-Newton, and 
is included in the XMM-Newton Serendipitous Source Catalog. 
The flux retrieved from the catalog is Fx = (7.75 ± 1.8) x 
10 _14 ergs _1 cm~ 2 . NGC 1860 is located close to the X-ray 
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Figure 2. Right panel: An ACIS-S image of the region of NGC 1818. produced from the event list binned by the factor of 8. The standard software detects 
NGC 1818 (within the square box), as well as other sources (encircled). Left panel: Enlarged area within the square box in the right panel with the adaptively 
smoothed image of NGC 1818. 



Table 2. Parameters of the massive star clusters 



Cluster Center (J2000.0) r c r c log Age Metallicity log M c \ log L^ s 

name a 8 [arcsec] [pc] [yr] [Fe/H] [Mq] [ergs ] 



R136 


05:38:42.5 


-69:06:03 


10 


NGC 2 100 


05:42:08.6 


-69:12:44 


5.02 


NGC 1818 


05:04:13.8 


-66:26:02 


10.10 


NGC 1850 


05:08:41.2 


-68:45:31 


10.48 


NGC 1860 


05:10:38.9 


-68:45:12 


9.00 


NGC 1831 


05:06:17.4 


-64:55:11 


18.28 


NGC 3603 


11:15:09.1 


-61:16:17 


29.5 


Arches 


17:45:50.4 


-28:49:20 


10 


Quintuplet 


17:46:14.8 


-28:49:35 


12 



2.43 


6.5 


~ -0.4 


4.40 


34.3 


1.22 


7.2 


-0.32 


4.31 


< 34.04 


2.45 


7.4 


0.0 to -0.4 


4.01 


34.32 


2.55 


7.5 


-0.12 


4.87 


< 32.70 


2.19 


8.3 


-0.52 


3.90 


34.61 


4.44 


8.5 


+0.01 


4.71 


< 33.48 


1.00 


6.0 


0.02 


4.0 


34.3 1 


0.39 


6.3 


0.02 


4.9 


34.2 2 


0.46 


6.6 


0.02 


3.8 


33.3 2 



Structural parame ters o f LMC clu sters are fromlMackev & Gilmorel 120031) 



1 from 

2 from 



Sung&Bessel<2004 and lMoffat et alj l2002t) 



Figeretal]|2002t) . Figer. McLean & Morris (1999) and lLaw & Yusef-ZadehU2004l) 



source 1AXGJ05 1054-6844 (see Fig.QJ. The latter is identified as 
a foreground late-type star in the literature. 

Based on the distance of the LMC , which we adopt as 50. 1 kpc 
throughout this paper l Mackev & Gilmorel2003l) . the X-ray fluxes 
corrected for the interstellar absorption are converted into luminosi- 
ties (Tabled. 

The most massive Galactic star clusters also have been ob- 
served in X-rays, and analysis of their X-ray properties is avail- 
able in the literature (Table|2j- Galactic clusters demonstrate the 
same trend as LMC clusters, with younger ones being more X- 
ray active. NGC 3603 (1 ± 1 Myr old) exhibits a wealth of stel- 
lar point sources, al ong with a detectable level of diffuse emission 
iMoffat et al.ll2002l) . The two most massive known in the Galaxy 
are the Arches (1-3 My r old) and the Quintuplet (3-5 Myr old) 
iFiger et alJll999tl2002l), located in the central 50 pc of the Milky 
Way. iLaw^^^jseL^adehl 120041) pointed out that the character of 
X-ray emission of the Arches and Quintuplet is different. The ratio 
of masses as well as of IR luminosities between Arches and Quin- 
tuplet is roughly 3:1. However, the ratio of X-ray flux (0.5-8 keV) 
is approximately 11:1; if only the point source emission is consid- 
ered, the ratio rises to 18:1. The Quintuplet cluster appears to have 



relatively more diffuse emission, the hottest gas being located in 
the cluster core. In contrast to the Arches' diffuse emission, the 
Quintuplet's diffuse emission is clearly thermal. 



3 MODELING OF X-RAY EMISSION FROM STAR 
CLUSTERS 

We consider an idealised star cluster with mass M c i = 1O 6 M . 
It is assumed that the stars are coeval, i.e. the cluster has under- 
gone an instantaneous burst of star formation. We assume that the 
stellar masses are distributed according to the standard (Salpeter) 
initial mass function (IMF), £(M) = M M~ 2 - 35 , with a lower 
mass cut-off at M m m = 1 Mq and an upper mass cut-off at 
A/m ax = 100 Mq. The radius of the cluster is such that the central 
density is about 10 5 Mq pc~ 3 . The calculations are done for solar 
(Z = 0.02) and LMC {Z = 0.008) metallicity. 

With the abov e assumptions, the output from the 
iLeitherer et alj 1 19991) stellar population synthesis code can 
be directly applied. 

The major simplifications we have made are the following: 
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(1) We effectively assume a cluster of single stars. The realistic 
account for binary stars in the cluster would affect the IMF and 
the age attributed to the cluster. We account for the binarity in a pa- 
rameterised way when modeling the cluster wind. We also consider 
binaries in estimates of the colle ctive stellar X-ray luminosi ty. (2) 
We neglect the cluster dynamics. iMackev & Gilmord 12003T) anal- 
ysed surface brightness profiles and structural parameters for 53 
clusters in the LMC. They concluded that the spread in the core 
radius increases significantly with increasing cluster age. For our 
model cluster the radius does not change during the evolution. Clus- 
ter mass is also assumed to be constant. (3) We assume that the 
mass-luminosity relation holds for the stars of all masses. We also 
assume that the stars evolve with constant mass till they reach WR 
phase. (4) We do not consider X-ray emission from degenerate stars 
and X-ray binaries. There is a broad discussion in the literature on 
this subject. (5) We exclude any possible sources of non-thermal 
X-rays. 

The time dependent X-ray luminosity of a cluster is 

Lx(i) = L x ° int (t) + Lx iff «. (1) 

The diffuse X-ray emission, L^ S , originates in the cluster wind. 
The collective stellar luminosity is 

L p x ° int (t)= / £(M) Lx{M,t) dM, (2) 
Ji/ nin 

where Lx(M,t) is the X-ray luminosity of a star with mass M 
at age t. In the next two sections we describe our modeling of the 
diffuse and the stellar component. 



4 DIFFUSE X-RAY EMISSION OF STAR CLUSTERS 
4.1 Wind from a model cluster 

The interstellar medium (ISM) in a star cluster receives a constant 
supply of gas, at a rate M*, due to the mass loss from massive stars 
and, later on, supernova explosions. In addition, material can be ab- 
lated from protostellar (or planetary) disks (e.g. lKroupa & BouvieJ 
120031) . Some quantities of primordial gas may also b e present in 
a young cluster. Following IStevens & Hartwelll 120031) . let us call 
this additional contribution M coo i. Then the total mass input to a 
star cluster is M = M* + M coo \. Stellar winds and SN explosions 
supply also kinetic energy to the cluster's ISM at a rate E. Some 
fraction of this kinetic energy input is thermalised, heating the gas 
to X-ray emitting temperatures. Both, M and E, vary significantly 
during evolution of the cluster. Being initially low for a very young 
cluster, stellar mass-loss and energy input are increasing while stars 
are evolving, and reach their maximum with the beginning of SN 
ignition. On the other hand, in an old cluster, where massive stars 
are absent, M and E steeply decline. 

Stellar populati on synthesis, such as Starburst99 
iLeitherer et alj Il999l) . are convenient tools to provide M 
and E that account for both, stellar winds and supernova explo- 
sions. We used here the standard output of Starburst99 for the 
appropriate metallicities to obtain the mass and energy yields. 

There is a shortcoming in Starburst99 for determining the 
mass and energy input in a dense star cluster. The model accounts 
only for single stars, while in fact the majority of high-mass stars 
in massive clusters are in binary or multiple systems. The mass- 
loss rate from a binary star is basically the sum of the mass- 
loss rates from two single stars. However, there can be a differ- 
ence in the rates of kinetic energy input between single and bi- 
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Figure 3. X-rays from a model cluster of 10 6 Mq . Upper panel: Luminos- 
ity of cluster wind versus cluster age for solar metallicity (solid line) and 
LMC metallicity (dashed line). The epochs when mass and energy input 
is dominated a) by stellar winds and b) by supernova explosions are indi- 
cated. The part of the curves between 5 and 35 Myr is artificially smoothed 
to correct for the computationa l discontin uity of Starburst99 models for SN 
rates (for details see lLeitherer et alll999l) . Lower panel: Temperature of the 
cluster wind as function of the cluster age. Solid and dashed lines are for the 
solar and LMC models, respectively. Note that after ss 40 Myr of cluster 
evolution, the temperature becomes too low for X-ray production. 

nary stars. In massive binary systems the stellar winds interact in a 
way, that a part of the total kinetic energy is dissipated, e.g. for 
the heating of the colliding winds zone, and radiated away. Al- 
though it is possible to speculate which fraction of kinetic energy 
is consumed due to the wind-wind collision, there are a number of 
other uncertainties regarding energy losses in a star cluster. A pa- 
ra meter a, accounting for all these uncertainties, was introduced 
in lCheval ier & Cle^gj jl98 5|) and late r called "thermalisation effi- 
ciencv" ir d^tevens&JHartwelll 120031) . Thus, the energy losses due 
to the wind-wind collisions in binary systems can be accounted for 
in a parameterised way, using factor a. 

To calculate the X-ray luminosity of cluster wind one 
must solve the hydrodynamical equations. IChevalier & Cleggl 
ll985l) have shown that a wind from a st ar formation region 
can be desc ribed as a self-similar process IZel'dovich & Raiseil 
11967. p.789l) . and obtained solutions applicable to mass loaded 
cluster winds. 

Let us consider the evolution of the broad band X-ray lumi- 
nosity of an idealised model cluster with mass 10 6 Mq and radius 
i?d = 2 pc. The resulting density is 3 x 1O 4 M0 pc -3 . We consider 
the X-ray luminosity of the whole cluster core, and do not infer its 
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Figure 4. Upper panel: Comparison of the observed diffuse X-ray luminosi- 
ties (diamonds) with individual cluster models. Crosses are the theoretical 
diffuse X-ray luminosity assuming that a = 1 and M coo \ = (see text). 
Circles are the model values assuming a = 0.1 and M coo i = 0.05JW*. 
Lower panel: Temperature of the cluster wind for a model with a = 1 and 
A^cool = (crosses) and a = 0.1 and M coo i = 0.05Af* (circles). 



radial profile. Therefore the lChevalier & CleggHl985t) solutions for 
density, p(r), and the pressure, P(r), are integrated over the cluster 
core radius. 

The frequency-integrated X-ray luminosity of the cluster wind 
at epoch t is 

Lf H {t) = A(T X ) EM(t), (3) 

where A(Tx) is the cooling function, and EM(t) is the emission 
measure of the cluster wind: 



EM(t) 



and 



Tx(i) 



dV, 



TTtHMi P(t) 
k p{t) ' 



(4) 



(5) 



with the Boltzman constant k and proton mass mu- The mean 
molecular weights per ion, p,j and per electron, p c , depend on the 
cluster age t due to chemical evolution. However, for simplicity 
we assume that pi and p c do not change significantly and consider 
them to be constant (1.3 and 1.14, respectively). 

From scaling laws of self-similar solutions, emission measure 
and temperature of the cluster wind within the mass and energy 
production region of radius 7? are 



EM(t) = 1.4 x 10" 



and 



4tt 



M 3 



/ii/i c m^ RaE 



k M 



(6) 



(7) 



Therefore, when M and E are known at time t, the temperature, 
density and velocity of the cluster wind can be obtained without 
solving partial differential equations. 

Assuming a — 1 and M coo \ = 0, we substitute the mass in- 
put, M(t), and kinetic energy input, E(t), provided by population 
synthesis into Eqs. OII^IEJ an d obtain the luminosity of the cluster 
wind versus time (Fig.[3j. 

As can be readily seen from Fig.|3| the level of diffuse emis- 
sion varies significantly during the cluster history. In a cluster 
younger than 2.4 Myr (log Age = 6.4), stellar population synthe- 
sis predicts neither stars with high mass-loss rates (such as LBVs 
or WR stars) nor supernova remnants. The influx of mass into the 
intracluster medium is therefore small, resulting in a relatively low 
density and emission measure. 

After the first ~ 2.4 Myr of stellar evolution the most massive 
stars pass trough the LBV stage, which is a very brief but neverthe- 
less important phase because of high mass loss. Subsequently they 
evolve to Wolf-Rayet stars. Due to the mass and energy input of 
these stars, the level of diffuse X-ray emission of the cluster wind 
increases by two orders of magnitude. WR stars finish their lives af- 
ter a few x 10 years by a SN explosion. 6.3 Myr (log Age = 6.8) 
after initial starburst the mass and energy input from SN explosions 
becomes dominant and maintains a nearly constant X-ray luminos- 
ity of the cluster wind for millions of years. 

It is not expected that stars with initial masses below 8 Mq 
undergo a SN explosion (Woosley, Heger & Weaver 2002). These 
stars do not have significant stellar winds either. Therefore, in a 
cluster older then ~ 40 Myr (log Age = 7.6) no stars with high 
mass loss are left. Consequently the mass and energy supply to the 
cluster wind is cut-off. The temperature drops sharply and cluster 
wind is not a source of diffuse X-ray emission any more. 



4.2 Comparison with observations 

The above theoretical considerations can now be confronted with 
the empirical data described in Section 3. First, we have plotted the 
observed diffuse X-ray luminosity of the clusters listed in Table|2|in 
the upper panel of Fig.[4](diamond symbols). Next, for each cluster 
we calculated the theoretical X-ray luminosity and the temperature 
of cluster wind, assuming a = 1, and M CO oi = 0. The crosses 
in Fig.|4] represent these theoretical values. As can be seen from 
the upper panel of Fig.|4|the observed diffuse X-ray luminosity is 
much higher than the theoretical one, the only exception being the 
Arches. 

Given the quality of the data, the information on the observed 
temperature is quite limited. lLaw & Yu sef-Zadeh 1 2004J) analysed 
the spectra of the diffuse gas in the Arches and in the Quintuplet. 
For the Arches, the best fit is a two-temperature thermal model with 
Ti = 1.4+°;;) x 10 7 KandT 2 = 7.2_ 6 . 6 x 10 7 K. For the Quin- 
tuplet the peak temperature of the diffuse gas in the clust er cen- 
ter is T = (3 ± 0.6) x 10 7 K. IStevens & Hartwell J2003I) report 
2 x 10 7 K as the temperature of diffuse gas in R 136. Our analy- 
sis of NGC 1818 suggests that the observed temperature is ranging 
from 4 x 10 6 to 6 x 10 7 K. (However, due to the quality of the 
spectra the question whether the Arches and NGC 1818 emission 
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is thermal remains open.) As can be seen from the lower panel in 
Fig-0the theoretical X-ray temperatures (crosses) are apparently 
too high compared to the temperature estimates deduced from the 
spectral fits. 

There are two aspects to be considered in order to explain 
these discrepancies. First, the luminosities of diffuse emission 
listed in Table|2| cannot be attributed unambiguously to the clus- 
ter wind alone. Second, regarding the model, we adopted a certain 
cluster age (listed in Table|2j. The estimates of the cluster ages, 
however, have a large uncertainty. E.g., the estimated ages are 1 - 
3 Myr for the Arches, 3 -5 Myr for the Quintuplet, and 1±1 Myr for 
NGC 3603. As can be seen from Fig.|3| the model predicts that the 
diffuse X-ray luminosity rises by two orders of magnitude between 
1 .6 and 3.2 Myr. Hence, the predicted X-ray luminosity is very sen- 
sitive to the cluster age. 

Second, the values represented by small crosses are obtained 
under the assumption that the mass-loss is only from stellar winds 
and SNe (M coo i = 0) and that all mechanical luminosity released 
in the cluster is spent to heat the cluster wind (a — 1). But, as 
can be seen from Eq. 0, a thermalisation efficiency a smaller then 
unity, and increased mass input will act to d ecrease the temper- 
ature of the cluster wind. As pointed out by IStevens & Hartwell 
120031) . this decrease of the temperature can be understood as a drop 
in the average energy available per particle. Since bremsstrahlung 
scales with the square root of temperature and using scaling rela- 
tion Eq. (7J, one notices that L x oc M 3 (ME)~ 1/2 . Therefore, for 
somewhat higher mass input and reduced energy input, the lumi- 
nosity is expected to rise. 

We did not attempt to infer Af coo i and a from the observations 
because of the many uncertainties in observed quantities. But it can 
be shown that it is possible to find such a combination of these pa- 
rameters that would make the theoretical values of X-ray luminos- 
ity and temperature resembling the observed ones. As an example, 
we plotted in Fig.|4| (large dots) Z/ X lfI and Tx that were obtained 
assuming plausible values A/ coo i = 0.05 M* and thermalisation 
efficiency a — 0.1. Compared to the models with A/ coo i = and 
a — 1 (crosses) the luminosity is higher while the temperature is 
lower. 

The theoretical values fit quite well to the observed ones for 
the R 136 and the Quintuplet. For the Arches, the model predicts 
higher diffuse X-ray luminosity than is actually observed. There 
are at least two explanations for this discrepancy. First, if the age of 
the Arches is actually smaller than assumed here (2 Myr), the pre- 
dicted X-ray luminosity will be reduced towards the observed one. 
Second, the diffuse luminosity of Arches may be higher than the 
one used here (1.6 x 10 34 erg s^ 1 ), dep ending on the interpretation 
of the observational data. For instance. [Stevens & Hartwell 1200 3l) 
used L x lff = 5 x 10 35 erg s _1 as the empirical value. 

Two LMC clusters, NGC 1818 and NGC 1860, show X-ray 
emission at a much higher level than expected from the cluster wind 
alone. To explain this, other sources of X-rays radiation should be 
invoked. From the image of NGC 1818 in Fig.|2|it is difficult to de- 
cide if the source is actually diffuse or pointlike. Also, the very low 
count-rate does not allow to extract any useful spectral information. 
However, one may expect the presence of a supernova remnant in a 
25 Myr old star cluster like NGC 1818. In fact our sample includes 
five LMC clusters older than 4 Myr. Statistically, at least one of 
them should exhibit a young supernova remnant (SNR). Applying 
scaling relations for the dynamical age of SNR and SN energy from 
Hughes, Hayashi & Koyama (1998) to NGC 1818, we estimate the 
possible SNR age of a fewxlO 3 years and SN's energy of about 
w 10 51 erg. The SNR X -ray emission further discussed in Sect. 6. 



The relatively high level of X-ray emission from the 200 Myr 
old cluster deserves special attention. The cluster was observed 9' 
off-axis. It is also 38" apart from a bright foreground star. As listed 
in the XMM-Newton Serendipitous Source Catalogue, NGC 1860 
has a detectable level of emission harder than 4.5 keV. As we will 
show in the subsequent sections, it is not expected that normal 
stars would be significant X-ray sources at the age of NGC 1818 
and NGC 1860. Therefore, X-rays in NGC 1860 can be possibly 
attributed to either a SNR or a degenerate star. 

Another interesting result is that the observed X-ray emission 
from NGC 3603 is much higher than can be accounted for wit h the 
cluster wind model. As was suggested bv lMoffat et all 120021) . the 
high luminosity of diffuse emission in this cluster may be attributed 
to unresolved low-mass stars. We will investigate the contribution 
of the low- and high-mass stars in the next sections. 



5 EVOLUTION OF THE X-RAY LUMINOSITY OF STARS 

If the X-ray luminosity of a star scales with its bolometric luminos- 
ity, and the latter scales with the stellar mass, one can estimate the 
total X-ray luminosity of stars with given mass. 

For simplicity we assume that stars evolve with roughly con- 
stant luminosity. In order to model the evolution of stellar X-ray 
luminosity of a cluster, we a) calculate the number of stars per 
mass bin from the adopted IMF; b) use a mass-luminosity relation 
to estimate the stellar bolometric luminosity; c) estimate the X-ray 
luminosity using a correlation between X-ray and bolometric lu- 
minosities; d) estimate the lifetime in hydrogen-burning stages for 
a star of given mass. Thus a correlation between X-ray luminosity 
and age is obtained. 

We have adopted the m ass-luminosity relation from 
iBinnev & Merrifieldil998. p.280h . 

f8l(^) 2 '" M>2OM 
Tf = { L78 (n|) 3 ' 8 2<M<20M Q (8) 
{ °- 75 (l^) 4 ' 8 M<2M Q , 

where Lboi is assumed to be equal to the bolometric luminosity on 
the main sequence. 

The origin of X-ray emission is different in low and high mass 
stars. Low-mass stars (M* < 3 Mq) emit X-rays as result of mag- 
netic activity caused by the dynamo mechanism and connected with 
the presence of outer convective layers. Massive stars, on the other 
hand, are thought to be X-ray emitters due to shock waves in their 
powerful (either colliding or freely expanding) winds driven by ra- 
diation pressure. Stars with masses between 3 and 10 Mq gener- 
ally do not emit X-rays, because outer convective zones are absent 
in stars more massive than 3Mq, and stellar winds are weak for 
stars less massive then 10 Mq. 

5.1 Low- and solar-mass stars 

To trace the evolution of X-ray emission from an ensemble of low- 
and solar-mass stars (0.1A4© < M* < 3Mq) we use the re- 
sults ob tained from the s t udy o f stars in the Orion Nebular Cluster 
(ONC). Fl accomio et alj 12003T) explored the relationship between 
X-ray activity and stellar mass and age for 696 well-characterized 
ONC members based on Chandra observations. They found a direct 
correlation between Lx and mass for stars less massive than 3 Mq . 
For stars with masses in the range from 3 to 10 Mq, Lx is small 
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Table 3. Depe ndence of X-ray lumin osity on age for low- and solar-mass 
stars (based on lFlaccomio et all2003l) 



M, [M Q ] 


0.5... 1 


1...3 


log (Age) [yrs] 


log(Lx) [erg/s] 


log(-Lx) [erg/s] 


5.5 - 6.0 


30.65 


31.50 


6.0 - 6.5 


30.0 


31.18 


6.5 - 7.5 


29.51 


29.73 




low-mass stars 



1 - 



6.0 6.2 6.4 6.6 6.8 

log( Time [yr] ) 

Figure 5. The collective X-ray luminosity of all cluster's low- and solar- 
mass stars (solid and dashed lines correspondingly) as a function of time 
(see text for details). For comparison the cluster wind luminosity is also 
shown (thin line). 



and attributed to the less massive companion in b inaries. For mas- 
sive O stars (M > 10 Mq), Lx rises again (see lFlaccomio et ail 
120031 Fig. 4). The authors attributed the drop in X-ray luminosity 
at ~ 3 Mq to the time when stars of the ONC's age (~ 1 Myr) 
become fully radiative and, therefore, the dynamo driven X-ray ac- 
tivity stops. The large number of ONC stars with very small masses 
(0.1 < M/Mq < 0.5) has been undetected in X-rays. 

In Table|3]we summarize the correlation between X-ray lumi- 
nosity and stellar age for stars with masses ( 0.5 < M/Mq < 3) 
observed in the ONC. lFlaccomio et alj 120031) discuss the interpre- 
tations of this correlations, including possible observational selec- 
tion effects. They conclude the existence of a real age spread in the 
ONC. 

Stelzer, Micela & Neuhauser (2004) observed the central re- 
gion of the ChamaeleonI South star-forming cloud with XMM- 
Newton. In general, they confirm the results obtained for the ONC 
and report that about 10% of the low-mass stars in their sample 
are X-ray flaring. Most of the variability can be attributed to the 
hard band (2.5-10 keV), while the emission in the soft band (0.2- 
2.5 keV) is more steady. Some of the variable stars seem to undergo 



irregular variability and flaring. A large number of stars can be de- 
tected only during short duration (~ 1 h) flares in the hard band. 

Let us trace the evolution of the total X-ray luminosity of low- 
and solar-mass stars in our 10 6 A/0 idealised model cluster 1 . The 
total number of stars with masses between 0.5 ... 1 Mq is 1.7x 10 , 
while 9.8 x 10 4 stars are in the 1 ... 3 Mq mass bin. 

Assuming that the luminosity of each star changes with age 
as specified in Table|3| and adding all stars in each mass bin, we 
obtain the collective X-ray luminosity of low- and solar-mass stars 
as function of age (Fig.|5}- Individual low- and solar-mass stars are 
relatively weak X-ray sources. However, they are numerous and 
difficult to resolve spatially. Figure[3] demonstrates that in clusters 
younger than 3 Myr the collective luminosity of low- and solar- 
mass stars is quite high. 

To understand the contribution of low-mass stars to the total 
X-ray emission of massive clusters let us imagine how the ONC 
would look if it were located in the LMC. Numerous X-ray dim 
low-mass stars fill the whole volume of the stellar cluster, while 
massive stars tend to concentr ate towards the center, Chandra's 
HRC observations analysed by iFlaccomio et alj 120031) covered a 
30' x 30' field of view, and were centered on the Trapezium clus- 
ter. The density of the Trapezium, ~ 5 x 10 3 stars pc -3 within the 
central 0. 1 pc, makes it one of the most dense clusterings known 
in our Galaxy. Ch andra observations resolved 111 point sources 
ISchulz et all2001T) in the Trapezium. While most of these sources 
are faint, twelve have an X-ray luminosity above 10 31 ergs -1 and 
four of those are identified with massive early-type stars. The ONC 
is located at 440 pc. At the distance of the LMC, the 807 stars from 
Chandra's HRC observations would be all within l'.'6 x l'.'6. The 
PSF of Chandra has a FWHM of 0'.'5. Hence it would be impossi- 
ble to resolve the low- and solar-mass stellar population, while the 
X-ray bright stars of Trapezium would appear as one point source 
embeded in diffuse emission. 

It was found in the previous section that the cluster wind emis- 
sion alone cannot explain the observed diffuse X-ray luminosity of 
NGC3603. This cluster is only about 1 Myr old, and one may ex- 
pect that stars less massive than 1. 5 Mb are still fully con vective 
and X-ray active (see references in lFlaccomio et al]|2003l) . Using 
the standard IMF and the average X-ray luminosity as listed in Ta- 
blcH we roughly estimate the collective X-ray luminosity of low- 
and solar-mass stars in NGC3603 as L^ w ~ 10 34 erg s~\ The 
sum of the cluster wind luminosity and the collective luminosity 
of low- and solar-mass stars matches the observed level of diffuse 
X-ray em ission in NGC 3603 ver y well, thus confirming the sug- 
gestion of lMoffat etall<2002h and lSung & Besselll 120041) . In older 
clusters, such as Quintuplet or R 136, the low-mass stars are less 
active, and massive stars become the main sources of stellar X-ray 
emission. 



5.2 O stars 

A correlation between X-ray and bolometric luminosity fo r Galac - 
tic O stars wa s sugg ested by ISeward & Chlebowskil ll982h . 
iBerghoefer et alj Il997t) compiled a catalog of optically bright 
OB stars (RASSOB) based on the Rosat All Sky Survey (RASS). 
They found that Lx ~ 10~ 7 Lboi for the en semble of O stars. 
From theoretical side. lOwocki & Cohenl 1 19991) explained this cor- 
relation for single stars invoking the wind-momentum luminosity 

For the estimate of the number of low-mass stars we assumed a lower 
mass cut-off of the IMF at 0. 1 Mq . 
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log(L bol /L Q ) log L bol / L 

Figure 6. Bolometric luminosity of O stars as function of their X-ray luminosity. Left panel: Spectroscopic binaries. Circles represent the X-ray lu- 
minosity for the Rosat passband (0.2-2.4 keV). Squares are for the XMM-Newton and Chandra data (0.2-10 keV). The line is the linear regression 
logLx = 1.05 log Lbol — 7.00. Right panel: Single runaway stars. The thick line is the linear regression logLx = 0.77 log L^ol — 5.94. The thin 
line is the linear regression obtained for the binaries. 



relation iKudritzki & Pulsll200Ch . We decided to verify this corre- 
lation selectively, distinguishing between binary and single stars 
and using the better data available now. 

Binary systems. We selected a sample of O+O spectroscopic binary 
systems with available Rosat pointing observations or observed by 
Chandra or XMM-Newton. In order to keep the data as uniform as 
possible, in case when observations by different instruments exist, 
the preference was given to Ros at (pointing observation s, in con- 
trast to the RASS data available to Berghoef er et alll997l) . The data 
are summarized in Table |4| If not indicated otherwise, the neutral 
hydrogen column density, distance and Lboi are taken from RAS- 
SOB. The count rates were retrieved from The Rosat Complete Re- 
sults Archive Sources for the PSPC catalog, and transformed into 
the X-ray flux assuming a Raymond-Smith thermal plasma with 
temperature kTx =0.6 keV. For the stars with resolved X-ray vari- 
ability, Lx was averaged over available observations. 

From linear regression analysis of our sample of spectroscopic 
binaries we obtain a good correlation log Lx ~ log Lboi — 7. This 
is rather puzzling because it is though that in binary stars a signifi- 
cant fraction of X-rays originates from a colliding winds zone. One 
may suggest that such a correlation is observed because Rosat's 
PSPC detector (0.2 - 2.4 keV) was sensitive primarily to the emis- 
sion from individual stellar winds, which have characteristic tem- 
peratures of about 0.6 keV, while plasma in the colliding wind zone 
might be heated to higher temperatures. If this were true, then the 
stars that were observed with XMM-Newton (0.2 - 10 keV) should 
not fit into the relation because this instrument also detects harder 
radiation. The fact that there is no significant scatter (see left panel 
of Fig.[6) between the binaries observed with XMM-Newton and 
Rosat makes us suggest that the bulk of observed X-rays originates 
in the individual stellar winds. 



|Pe Becker etal J J20041 concluded from their study of the col- 
liding O+O binary HD 159176 that the bulk of X-rays comes from 
a plasma with kTx = 0.6 keV, as expected from the hydrody- 
namic models of single stars, and that the line width is consistent 
with the stellar wind velocity. However, the steady-state colliding 
wind model systematically predicts too high X-ray luminosities. 
Nevertheless, the observed "hard tail" in the X-ray spectrum of 
HD 159176 most l ikely originates from the colliding wind zone. 
iMiller et alj 120021) discussed extensively the origin of X-rays in 
the eclipsing system 8 Oh. They as well concluded that the bulk of 
X-rays is originating in the O star wind. 

Single runaway stars. If our hypothesis about stellar wind emis- 
sion dominance over colliding wind zone emission in O binaries is 
correct, then there should be a similar Lx oc Lboi correlation for 
single stars. There is no way to select a sample of O stars which are 
definitely single. However, one might select a sample of runaway 
stars. The fraction of single stars among all runaway stars is very 
high (~ 75%). We hav e chosen all single runaway stars from the 
O Star Speckle Survey iMason et alll998l) (see Table and per- 
formed a linear regression analysis for this sample. There are only 
seven single runaway stars with available Rosat pointing, therefore 
the result is not statistically significant. Nevertheless, the X-ray lu- 
minosity of single stars appears to correlate with Lboi in a similar 
way as found for the O star binaries (Fig. [6). 

Hence we may apply the correlation log Lx ~ log Lboi — 7 
for O stars, disregarding whether they are single or binaries. 

5.3 LBV stars 

A massive star with initial mass above 50 Mq that began its life on 
the MS as an O type star may become very unstable near the end 
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Table 4. Pointing X-ray observations of O-type spectroscopic binaries and single runaway stars 



Name Exp. Time Count Rate log Nu d log Lx log -Lbol Detector 

[ksec] [10 2 cts _1 ] [cm -2 ] [pc] [ergs -1 ] [ergs -1 ] 



Spectroscopic binary stars 



HD 57060 


3.1 


3.55±0.4 


20.81 


1514 


32.15 


39.44 


PSPC 


HD 57061 


77.7 


10.93±0.1 


20.74 


1514 


32.62 


39.46 


PSPC 


HD 93205 


1.4 


5.90 


21.33 


2630 


33.04 


39.50 


PSPC 


HD 93206 


1.2 


6.59±0 


21.34 


2512 


33.06 


39.46 


PSPC 


HD 93403 


37.8 


2.10 


21.57 


3200 


33.57 


39.70 


EPIC 


HD 97484 


4.3 


1.72±0.2 


20.44 


3200 


32.41 


39.28 


PSPC 


HD 135240 


2.5 


7.30±0.6 


21.18 


615 


33.59 


38.77 


PSPC 


HD 149404 


2.8 


8.26±0.6 


21.40 


1380 


32.68 


39.70 


PSPC 


HD 152218 


7.6 


l.llztO.l 


21.34 


1905 


32.05 


38.78 


PSPC 


HD 152219 


7.6 


2.21±0.2 


20.56 


1900 


31.98 


38.76 


PSPC 


HD 152248 


180.0 






1757 


33.00 


39.81 


EPIC 


HD 152590 


8.9 


0.91+0.1 


20.57 


1900 


31.85 


38.65 


PSPC 


HD 159176 


20.4 


130.30±0.9 


21.33 


1500 


33.17 


39.17 


EPIC 


HD 165052 


3.5 


10.41±0.6 


21.36 


1585 


32.88 


39.10 


PSPC 


HD 167771 


5.9 


4.2±0 


21.40 


1585 


32.54 


38.97 


PSPC 


HD 203064 


6.5 


5.1±0 


20.98 


794 


31.79 


39.09 


PSPC 


HD 206267 


3.1 


10.4±0.6 


21.49 


780 


32.36 


39.25 


PSPC 


<50riA 


2.3 


176.5±2.8 


20.18 


501 


32.77 


39.34 


ACIS-S 


tOri 


0.7 


158.9±4.9 


20.30 


501 


32.68 


39.18 


PSPC 


LZCep 


5.0 


5.64±0.3 


21.33 


583 


31.72 


38.66 


PSPC 


V729Cyg 


3.4 


12.20±0.6 


19.94 


1800 


32.96 


39.7 


PSPC 


Plaskett 


3.7 


15.60±0.7 


21.18 


1514 


32.91 


39.35 


PSPC 


15Mon 


0.9 


3.00±1.8 


20.35 


692 


32.31 


38.89 


PSPC 



Runaway single stars 



£ Per 


1.0 


20.00±1.0 


21.06 


398 


31.68 


39.47 


PSPC 


a Cam 


3.8 


6.59±0.4 


21.08 


1010 


32.10 


39.42 


PSPC 


/iCol 


3.7 


15.14±0.7 


20.04 


669 


31.89 


38.45 


PSPC 


AE Aur 


2.4 


4.37±0.4 


21.30 


321 


31.02 


38.16 


PSPC 


COph 


2.0 


39.00±2.0 


20.78 


154 


31.15 


38.46 


PSPC 


ACep 


4.2 


6.50±0.6 


21.20 


832 


31.92 


39.42 


PSPC 


CPup 


55.7 


126.90±1.6 


20.00 


437 


32.43 


39.49 


PSPC 



count rates are from "Complete Results Archive Sources for the PSPC" catalog 
unless otherwise indicated 



HD93403: L x fromlRauw et al. 


2002! 




HD152248: L x from 
HD159176: L x from 
V729Cyg: L bol from 


Sana et al 


2004 




De Becker et al. 


2004 


iRauwetalJl 19991) 



of core H-burning. Such stars are called luminous blue variables 
(LBVs) or S Dor variables. They populate the empirical upper lumi- 
nosity boundary in the H-R diagram and evolve directly to the WR 
stage via high mass loss episodes. Less luminous LBVs have initial 
masses less than 40 Mq and have probably already been through 
the red supergiant stage faumpnrevsll20(m Dozens of LBVs are 
identi fied in the Milky Way and the Local Group fvan Gendered 
1200 ll) . Many LBVs are found in the vicinity of a young star cluster. 

Some of the LBVs are luminous X-ray sources. iNaze et ail 

120021) deduce from the Chandra observation of NGC 346 in the 
SMC that the X-ray luminosity of the binary LBV star HD 5980 is 
1.7 x 10 34 erg s _1 . This in fact exceeds the total luminosity of the 
NGC 346 cluster (1.5 x 10 34 ergs -1 ). The famous Galactic LBV 
77 Car has a similar luminosity Lx = 1.2 x 10 35 erg s _1 (Leuteneg- 



ger, Kahn & Ramsay 2003). IPittard & Corcorai] <2002t fitted the 
high resolution X-ray spectrum of 77 Car with hydrodynamic mod- 
els of colliding winds. The X-ray emission of i) Car is similar to 
WR 140 and 7 Vel: the emission is hard and variable, with the min- 
imum near the phase when the star with denser wind is presumably 
in front. 

However, not all LBVs are bright in X-rays. The prototypical 
LBV star PCyg was observed for 97 ksec with R osat. Using in- 
terstel lar reddening and bolometric luminosity from Ivan Gendereil 
1200 ll) we estimate the X-ray luminosity of PCyg as Lx = 
2 x 10 31 ergs _1 resulting in logLx/Lboi = —8. Another ex- 
ample is the blue supergiant Sher25 asso ciated with NGC 3603 
iBrandner et alll997t) . lMoffat et ai"ll2002l) report that this star has 
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no significant X-ray detection and attribute this to the relatively low 
Lboi combined with a low wind velocity. 

5.4 WR stars 

During its evolution a massive star looses a significant amount 
of mass via a stellar wind, LBV outbursts or mass trans- 
fer in a close binary system, revealing first the CNO-burning 
products at its surface (WN phase), and subsequently the He- 
burning products (WC phase). The mass-loss from WR stars 
(M W r ~ 1CT 5 - _6 M Q yr" 1 . iHamann & Koesterkdfl993) is at 
least one order of m agnitude higher than in O-type stars (Mo ~ 
10" 6 - " 7 M Q vr~ 1 . lRepolust et alfeool) . The winds of WR stars 
are enriched by nuclear evolution products. An understanding of 
the mechanisms by which radiative pressure drives W R stellar 
winds is currently emerging jGrafener & Hamanrl2004l) . 

X-ray emission from single WR stars is largely enigmatic. 
There are no theoretical models of WR winds that describe the gen- 
eration of X-rays, although the expectation is that the same mecha- 
nism as in O winds may operate here. T his is supported by a stud y 
of instability growth rates in WR winds iGavlev & Owockil 19951) . 
Meanwhile, the increasing number of observations of WR stars re- 
veal a complex picture. Pollock, Haberl & Corcoran (1995) anal- 
ysed the RASS observations of WR stars, and concluded that there 
is no correlation between bolometric and X-ray luminosity. 

Similar to the study of O-type stars we compile the most 
sensitive up-to-date X-ray observations of WR stars (Table|5j. We 
restric t ourselves to WN-type stars. As shown in lOskinova et alj 
120031) . WC-type stars were not detected in X-rays, except for a 
few binaries. Figure^ shows the bolometric and X-ray luminosi- 
ties of our sample of Galactic WN stars. We first discuss the X-ray 
emission from putatively single WN stars, and then from binary 
stars. 

Single WN stars. Almost all bright O stars are found to emit X- 
rays. Out of the 1 1 stars in the RASSOB sample which were not 
detected in the Rosat All Sky Survey, nine have been detected since. 
Only HD 68450 and HD 154368 still remain undetected, probably 
because being observed only off-axis with short RASS exposures. 

In contrast, the X-ray luminosity of WN stars shows a much 
larger scatter. While a couple of WN stars have significant X-ray 
fluxes, some remain below the detection limit even after long expo- 
sures. Correcting for their distance and reddening, which are typ- 
ically both larger than for the RASSOB stars, these WN stars still 
must be intrinsically X-ray dim. 

Most intriguing are the XMM-Newton observations of WR61 
and WR40, and the Rosat pointing observation of WR 16. These 
stars were not detected, and upper limits on log (Lx/Lboi) are - 
8.1, -9.2 and -8.9 for WR61, WR40 and WR 16, respectively. 

Note that WR40 and WR 16 belong to the spectral subtype 
WN8. Remarkably, no WR star of this spectral type was ever de- 
tected in X-rays (see Table|5}. Chand ra HRC-I direct im aging of 
the binary WR147 (WN8(h)+B0.5V) <Pittard et alj|2002h showed 
X-rays which are certainly not co-spatial with the WN8 compo- 
nent, but associated either with the companion star or the colliding 
wind zone. There are indications that stellar winds from WN8 stars 
are basically different from other WN-type stars (Grafener, priv. 
communication) . 

As can be seen from Fig.0 the scatter of X-ray luminosity for 
the detected WN stars is extremely high, reaching up to three orders 
of magnitude. Apart from the WN8 case discussed above, there 
is no correlation with the WN subtype. The outstandingly X-ray 
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Figure 7. X-ray versus bolometric luminosity of WR stars observed by 
Chandra, XMM-Newton, or with Rosat PSPC pointings. Diamonds denote 
putatively single WN stars and circles denote WN spectroscopic binaries. 
Numbers are for the stars from the Vltth Catalog of WR stars. Upper limits 
for non-detections are shown by the arrows. The data are from Tablel5l The 
thin straight line represents the correlation between Lx and Lboi as found 
for the O-stars SBs. 

bright star WR 25 was carefully examined in lRaassen et alj 120031) . 
who proposed that WR 25 might be a binary in order to explain its 
rather hard X-ray spectrum. 

Binary WN stars. We have examined the available X-ray obser- 
vations of WR binary stars. One half of 22 spectroscopic binary 
WN stars listed in the Vllth Catalog of WR stars was not detected 
by RASS (the detection rate is even lower for WC-type binaries). 
Among those which were detected, we select the stars with the most 
sensitive observations (Table|5}. As can be seen from Fig.0 the 
X-ray luminosity is proportional to the total bolometric luminos- 
ity of the system. Moreover, the correlation is the same as found 
for O+O binaries, log Lx ~ log Lboi — 7. The only exception is 
WR 141, which was in Chandra's field of view but inconveniently 
located at a CCD edge. The Chandra source detection software 
recognises WR141. The count rate is 1.4 x 10~ s cts -1 , yield- 
ing log(Lx/Lboi) ~ —8.7, which is the lowest value among all 
detected massive binaries. Therefore we consider this count rate 
being uncertain and exclude the star from our statistical sample to 
probe the Lx versus Lboi relation. 

One should bear in mind that in general the O star is the more 
luminous component in a WR+O binary. Similar relations between 
Lx and Lboi in WN+O and O+O binaries can be expected if the 
major fraction of observed soft X-rays originates in the individ- 
ual stellar winds rather than in the colliding winds z one. It seems 
that th is preposition can be confirmed observationallv. lMaeda et all 
J 19991) analysed ASCA observation of WR139 (V444Cyg). They 
attribute the soft-component emission at kTi ~ 0.6 keV and non- 
variable luminosity to the individual 06 and WN5 components of 
the system. The hard component (kT2 ~ 2 keV) is phase variable 
and is caused by a colliding wind shock. Our preliminary analy- 
sis of XMM-Newton observations of WR 22 gives a similar result. 
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Table 5. Sample of X-ray observations of WN-type stars 



WR 


Sp. type 


log L bol 


lo g 7V H 


d 


Exp. 


logL x 


Detector 






[ergs -1 ] 


[cm" 2 ] 


[kpc] 


[ks] 


[ergs -1 ] 




Putatively single WN stars 


2 


WN2w 


39.24 


21.30 


2.51 


9.2 


31.9 


PSPC 


3 


WN3 


39.19 


20.96 


4.30 


8.4 


32.4 


EPIC 


152 


WN3w 


38.84 


21.26 


2.75 


4.7 


31.5 


PSPC 


6 


WN4s 


38.99 


20.78 


0.97 


9.5 


32.5 


EPIC 


1 


WN4s 


38.89 


21.48 


1.82 


7.6 


32.3 


EPIC 


7 


WN4 


38.89 


21.33 


5.75 


7.0 


32.3 


PSPC 


18 


WN5s 


39.29 


21.47 


4.57 


7.6 


32.4 


PSPC 


42d 


WN5 


39.29 


21.85 


4.09 


49.3 


32.6 


ACIS-I 


44a 


WN5 


39.29 


21.91 


3.24 


57.5 


31.2 


ACIS-S 


110 


WN5-6 


39.49 


22.02 


1.28 


20.5 


32.6 


EPIC 


20b 


WN6 


39.18 


21.94 


2.27 


36.2 


32.6 


ACIS-I 


136 


WN6 


39.69 


21.36 


1.26 


93.9 


30.9 


ACIS-S 


78 


WN7 


39.78 


21.28 


1.99 


9.9 


31.1 


PSPC 


25 


WN7a 


40.38 


21.28 


3.24 


109.0 


33.9 


EPIC 


79a 


WN9 


39.59 


21.28 


1.99 


9.4 


31.9 


PSPC 


Binary WN-type stars 


46 


WN3+OB? 


39.09 


21.08 


4.07 


70.3 


32.5 


EPIC 


133 


WN4+09I 


39.36 


21.15 


2.14 


78 


33.0 


GIS+SIS 


138 


WN5+B? 


38.88 


21.38 


1.26 


5.43 


32.2 


PSPC 


139 


WN5+06III-V 


38.9 W r + 39.1o 


21.50 


1.90 


19.0 


32.5 


PSPC 


141 


WN5+05V-HI 


39.97 


21.66 


1.26 


19 


30.6 


ACIS-S 


47 


WN6+05V 


39.60 


21.65 


3.80 


151 


33.05 


EPIC 


155 


WN6+09II 


39.08 W R + 39.29o 


21.40 


2.75 


90 


32.64 


GIS+SIS 


22 


WN7h+09III-V 


40.14 


21.02 


3.24 


82.7 


33.17 


EPIC 


20a 


WN6+WN6 


39.99+39.99 


21.89 


5.75 


32.6 


33.46 


ACIS-I 


147 


WN8(h)+B0.5V 


39.68 


22.12 


0.65 


4.9 


32.91 


HRS-I 


Significant non-detections of WN stars 


61 


WN5w 


38.79 


21.41 


6.79 


52.6 


< 30.7 


EPIC 


157 


WN5 


39.49 


21.59 


3.39 


5.8 


< 31.6 


PSPC 


40 


WN8 


39.29 


21.24 


2.26 


16.5 


< 30.1 


EPIC 


16 


WN8h 


38.85 


21.36 


2.37 


7.6 


< 30.0 


PSPC 


124 


WN8h 


38.88 


21.70 


3.36 


4 


< 32.3 


HRI 


148 


WN8h+B3IV/BH 


39.59 


21.62 


8.28 


9.9 


< 32.2 


IPC 



X-ray luminosities are based on: 
PSPC count rates from "ROSAT Complete Results Archive Sources for the PSPC" catalog unless otherwise indicated 
EPIC count rates obtained from the analysis of XMM-Newton archive pipeline products unless otherwise indicated 
ACIS, HRS count rates obtained from the analysis of Chandra archive pipeline products unless otherwise indicated 

WR 1 : Lx from Igna ce, Oskinova & Brown ( 2003) 

WR 3 : Distance fromlMarchenko et alj <2004l) 

WR6: L x from lSkinner et all l2002bT 

WR25: L x from lRaassen et alj<2003l) 

WR20a: Spectral type from lRauw et alj<2004 . L bo i is assigned according to the spectral type and mass 

WR 42d, 44a, 79a: L br ,i is assigned a ccording to the spectral type 

WR46: L bo i from lVeen etai]l2002h 

WR79a: Marginal detection at 12min of f-axe 

WR 110: L x from lSkinner et al J 12002 J) 

WR 139: L bol from lMarchenko et alJi 199^1 

WR 141 : Located at a CCD edge 

WR 147: Distance and extinction from lPittard et alJ i20O%) 

WR 133, 155: X-ray luminosity is roughly estimated from ASCA data 

WR7, 16, 18, 124, 136, 141, 148, 157: distance, L b _ v , and L bol are from lHamann & Koesterkehl998h 
WR22, 25, 141, 147, 155: Stellar parameters from Hamann et al. (in prep.) 

WR 133, 138: L bol from Harries, Hillier & Howarth (1998). For WR 133 we add a typical luminosity of an OI star log L bol /L Q = 5.7 
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Overall, given the small number of available observations, it ap- 
pears that WN+O binaries are quite similar to O+O binaries with 
respect to a correlation between Lx and Lboi- 

Summarizing, we conclude that the major fraction of X-ray 
photons in colliding wind binaries is emitted by the individual stel- 
lar winds of the binary components. The individual stellar winds 
typically have thermal spectra with characteristic temperatures of 
fcTx ~ 0.6 keV. The colliding wind zone may manifest itself by 
the presence of an additional, somewhat harder, either thermal or 
non-thermal component. 

Binary WC stars. Only six WC spectroscopic binaries have been 
observed in X-rays with exposures longer than that of the RASS. 
ASCA observations of WR 132 and WR 1 13 are still waiting to be 
carefully examined. WR 1 1 (7 Vel, WC8+07.5 III) was extensivel y 
observed with Chandra and XMM-Newton ("e.g. lSchild et alj2004l) . 
From an analysis o f the spectra taken in different orbital phases, 
ISchild et ail 1 2004) conclude that the total X-ray emission of the 
system is dominated by the material located in the wind collision 
zone. The X-ray emission from the wind of the O-type companion 
is one order of magnitude smaller, and from the WC-type com- 
panion is negligible. For this system logLx/Lboi = —5.9, vary- 
ing by w 0.5 dex in dependence on orbital pha se. ASCA observa- 
tions o f WR 140 (WC7+04-5V) are reported bv lZhekov & Skinned 
l200Ch . This star has also been monitored by RXTE and Chan- 
dra. There is no detailed stellar atmosphere analysis available 
for WR 140. We roughly estimate the bolometric luminosity of 
the system from the spectral types of the components, adopting 
Lboi(WC7) w 2 x 10 5 L o and L bo i(06V) w 1.3 x 10 5 L Q . 
Then the ratio of X-ray and bolometric luminosity of the system 
is estimated as log L^/L hol « -5.8. WR48 (6 Mus, WC6+O6V) 
was observed by a Rosat PSPC pointing. Assigning the same bolo- 
metric luminosity as for WR140 we find logix/iboi ~ —6.1. 
WR79 (WC7+05V) was observed by XMM-Newton. We have re- 
trieved and analysed the archival data for this binary. The spec- 
tral types of the companions are similar to WR 140 and WR48, 
therefore we use the same bolometric luminosity for an order of 
magnitude estimate, and obtain logLx/Lboi ~ —7.4. WR79 is 
located in the young star cluster NG C 6231, where a large number 
of massive binaries is present iGarcia & Mermilliodll200ol) . Inter- 
estingly, WR79 is the brightest X-ray source in the cluster in the 
4.5-12.0 keV band. The spectral energy distribution of WR79 is 
quite hard. Perhaps, by analogy with 7 Vel, this hardness is due to 
the strong absorption that soft X-ray photons suffer when passing 
through the opaque WC wind. 

From our brief review of WC binary systems it appears that 
the X-ray emission from the colliding wind zone dominates over 
the emission of the individual stellar winds. Our tentative analysis 
does not reveal any Lx-Lboi correlation for WC binaries. 

Population synthesis predicts that most of the time during 
cluster evolution, the number of WC stars is much smaller than 
of WN stars. However in clusters that are about 3-4 Myr old, the 
number of WC stars c an be comparable to the number of WN stars 
iLeitherer etalll999h . 

Regarding to the X-ray emission from WR stars we may con- 
clude that: (a) the emission from binary systems does not signif- 
icantly differ between O+O and O+WN systems; (b) the level of 
emission from putatively single WN stars does not correlate with 
the bolometric luminosity of the stars. Some single WR stars are 
quite luminous X-ray sources, while others appear as non-emitters; 
(c) Although some of the WC binary systems can be relatively X- 
ray bright, this is not universal for all WC binaries. No X-ray emis- 
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Figure 8. Lifetime in the H-burning stage, tms, as function of the initial 
mass, M zamB , for stars with solar metallicity. The results from lHirschi et alj 
|2004|) are for models with rotation (f ro t = 300 km s —1 , d ashed line) or 
without rotation (dotted). Models from lSchaller et at] ll992t) . also without 
rotation, extend to higher masses (solid line). 



sion yet detected from single WC stars. Hence, it seems that Lx 
decreases for a significant fraction of the massive stars when they 
evolve to the WN stage, and drops even lower in the WC stage. 

5.5 Evolution of X-ray luminosity of massive stars 

The predicted stellar lifetimes depend on the stellar evolu- 
tion model. Recent Geneva models include the effects of rota- 
tion (Hirschi, Meynet & Maeder A., 2004). The effects of bi- 
narity on the evolution are in cluded in the models by, e.g., 
I Van Bever & Vanbeverenl J2003). Both effects extend the lifetime 
on the main sequence. Therefore one attributes a younger age to 
a cluster when applying these models. Figure|8| shows the life- 
ti me on the H-b urning stage f or stars of different m ass as given 
bv lHirschi et"ai] <2004l) and by ISchaller elaH <1992l) . We use the 
latter lifetimes in this paper. We do not explicitly include LBVs 
because of their fast evolution (~ 40 000 yrs), and because little 
is known about their X-ray emission. However, an LBV star may 
be the do minant X-ray sou rce in a stellar cluster, as observed in 
NGC 346 iNaze et all20o3) . 

Figure|5| shows the evolution of the collective X-ray luminos- 
ity of massive stars in our model cluster. Considering Fig.|9|we shall 
notice an interesting interplay between the number of massive stars 
and the level of diffuse emission from the cluster wind. The cluster 
wind is fed my massive stars. If there are iV„ mass-loosing stars in 
the cluster, the cluster wind luminosity L x lfI oc M 3 cx (see 
Eqs.[3][6j while the collective stellar luminosity L x oc AT*. There- 
fore, the ratio of cluster wind to stellar emission is smaller for less 
massive clusters. 

Another point to address is the relation between the X-ray 
luminosity of low- versus high-mass stars. Our model predicts 
that in very young clusters the collective luminosity of low-mass 
stars should exceed the collective luminosity of the high-mass stars 
(compare Fig.|5|and Fig.[5). Nevertheless, the massive stars are de- 
tected as point sources in, e.g., NGC 3603, and are not out-shone by 
low-mass stars. The massive stars tend to concentrate in the clus- 



14 


L. M. Oskinova 


4 


1 i 1 i I i I i I i 

diffuse emission 


£ 3 

X 

>-J 

"Hi 
_o 


/ 

massive stars / 


2 






i.i.i.i,!, 


6.0 6.2 6.4 6.6 6.8 
log( Time [yr] ) 

Figure 9. Evolution of the collective X-ray luminosity of massive stars in 



a cluster M c \ = 1O 6 M0 as function of age. The X-ray luminosity of the 
ensemble of massive stars strongly declines with age. 

ter core, while low-mass stars fill the whole volume of the clus- 
ter. Therefore the surface X-ray brightness due to the low-mass 
stars is low, while the X-ray bright massive stars can be easily 
detected above the background. This situation lasts till the most 
massive stars enter the WR phase, and start to actively feed the 
cluster wind. After 2 .. 3Myr of cluster evolution, the emissivity of 
the cluster wind exceeds the emissivity of individual stellar X-ray 
sources. This is exactly the situation observed in the LMC clus- 
ters, and shown in Fig.Q the massive stars are seen in the youngest 
cluster, but sink into oblivion in the older clusters. 
As seen in Fig.|51 the x 

-ray luminosity declines steeply after 
w 2 Myr of cluster evolution, when the most massive stars end their 
lives by a core collapse and a supernova explosion. Black holes 
are formed from stars with initial masses above ~ 21 Mp,, while 
neutr on stars result from stars with lower masses IWooslev et alj 
120021) . From the analysis of the observations it appears that ~ 2 — 
5% of all OB stars must produce high mass X-ray binaries (HMXB) 
(consisting of a compact object and a massive s tar). This estimate 
is con sistent with population synthesis studies iHelfand & Moranl 
l200ll) . The lifetime of a HMXB is limited by the nuclear time-scale 
of the optical companion (~ 20 Myr). 

As can be seen from Tabled NGC2100, NGC 1818 and 
NGC 1850 are old enough to produce compact companions in 
HMXBs, and are young enough for most non-degenerate compan- 
ions still to be present. Adopting the quoted estimate of ~ 2 — 5% 
HMXB from all OB stars, one should expect 1 - 3 HMXBs in each 
of these three clusters. However, none of them shows X-ray point 
sources which are bright enough that they could be attributed to X- 
ray binaries. This apparent contradiction may indicate a problem 
with the predicted HXMB formation probability. 



6 X-RAY EMISSION FROM SUPERNOVA REMNANTS 

There are distinct phases in the dynamics of supernova rem- 
nants. Initially, the energy liberated in the collapse is deposited 
in the stellar envelope. The envelope is heated to a high temper- 



ature and ejected with high velocity. The expansion is uniform, 
with v oc r till the mass of swept-up material becomes sig- 
nificant. The highly supersonic expansion shock-compresses the 
ISM. The shoc ked heated ISM is hea ted and radiates by thermal 
bremsstrahlung. lBregman et al. U200H) analysed the ROSAT obser- 
vations of the complete sample of nearby supernovae that occurred 
between 1985.5 and 1994.3. They concluded, that the probability of 
an individual SN to have a luminosity higher then 2 x 10 39 erg s _1 
is less then 12%. The probability of a supernova to be more lu- 
minous than 6 x 10 3 8 erg s -1 is in the range 9% — 51%. From 
iLeitherer et al] j 19991) the SN rate for a solar metallicity cluster is 
about 10~ 3 SNyr" 1 for a 10 6 M W cluster older than ~ 3 Myr. Ap- 
plying [Bregmaj^Uii] l2003h results, the rate of X-ray bright su- 
pernovae in such a cluster can be roughly estimated as ~ 10 -4 per 
year. When an X-ray bright SN occurs, it is likely to be the brightest 
source of X-rays in its cluster, but only for rather short time. The 
X-ray brightness of SNe declines sharp with time as t -1 or t~ 2 . 

After a few x 10 2 years, the swept-up mass becomes greater 
then the ejected mass, and the expansion is described by the adia- 
batic blast-wave (Sedov-Taylor) solution. The radius of the remnant 
is given by 

r s = 0.31(E/n o ) 1/5 t 2/5 [pc], (9) 

where E is the SN explosion energy in units of 10 erg, no is the 
number density of the preshock ISM, and t is the SNR age in years 
and it is assumed that the material has solar composition. The post- 
shock temperature is 

kT s = 1.8 x 10 5 [keV]. (10) 

The frequency-integrated X-ray luminosity of the remnant is 
L X NR = EM Ax- The emissivity, Ax, depends on T s , electron 
temperature (T c ), abundances, and the ionisation state of the gas. 
For strong shock conditions and ideal gas behavior, the postshock- 
to-preshock density ratio is constant, n/no — 4. With the depen- 
dence of the remnant radius on time in Eq. |9j> the emission mea- 
sure of the remnant scales as EM oc rigt 6 ^ 5 . 

Borkowski, Lyerly & Reynolds (2001) presented improved 
calculations of the X-ray spectra for SNRs in the Sedov-Taylor 
phase. Their model is implemented in the xspec software package. 
The input parameters of this model are the shock temperature T s , 
electron temperature T c and the ionisation timescale, determined 
by the product of post-shock electron density and the remnant's 
age. Using the cluster wind model to obtain pre-shock density and 
T s (see Eq.[l0}, we could, in principle, model SNR X-ray emission 
for a cluster of given age, and trace how the environment in the 
cluster affects the SNR spectra. However, XSPEC models are not 
implemented over a wide enoug h range of parame ters. 

Using formulae from lBorkowski et alj 1200 lh and our cluster 
model, we calculated the time it will take for a SNR in Sedov- 
Tayler phase to cool down below 0.2 keV. It was assumed that 
T c = 0.5 T s and E = 10 51 erg s _1 . The result is shown in Fig.fTol 
As can be seen, the cooling time is slowly increasing for older clus- 
ters. However, even taking into account this small effect, in a dense 
cluster, the SNR remains a significant X-ray source, with unab- 
sorbed L| NR « 10 35 " 36 ergs -1 , for only a few thousand years. 
During this time the L X NR is comparable to the level of the cluster 
wind luminosity (see Fig.[3) and the remnant radius is compared to 
the cluster core radius. However, since the cluster wind temperature 
is high and roughly constant (see lower panel in Fig.|5J, the SNR 
spectrum is likely to be softer than the cluster wind spectrum. The 
supernova rate between 5 and 35 Myr in a 10 6 A/q cluster is about 
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Figure 10. Time for the shock temperature kT B of a SNR in Sedov-Tayler 
phase to drop below 0.2 keV, as function of the model cluster age. The clus- 
ter has solar metallicity. 

10 _3 yr _1 lLeithereretailll999h . therefore we may expect an en- 
hanced level of diffuse emission for older clusters, as it is observed 
in the LMC (see Fig.|3J. 

During the further phase in the dynamics of a SNR, as the 
remnant continues to expand, the cooling by optical emission lines 
become important and X-ray emission drops. Later the SNR expan- 
sion becomes subsonic, and the remnant disperses into the ISM. 



7 CONCLUSIONS 

We model both the thermal X-ray emission originating in a cluster 
wind and the evolution of X-ray emission from a stellar population 
and compare the model with observations. The conclusions are: 

(1) Effects of stellar evolution in massive clusters can explain the 
observed differences in the level of diffuse emission and numbers 
of point sources in clusters of different age. 

(2) The observed level of diffuse X-ray emission in young massive 
star clusters is in accordance with the theory of mass-loaded cluster 
winds. We retrieved and analysed X-ray observations of six promi- 
nent LMC clusters. The striking differences in the level of X-rays in 
these clusters can be explained by the different evolutionary stages 
of the clusters. The same holds for massive Galactic star clusters: 
NGC 3603, the Arches and the Quintuplet. 

(3) The X-ray luminosity of the cluster wind is low for the youngest 
clusters, such as NGC 3603, where the majority of massive stars 
is still on the main sequence and, therefore, have moderate stellar 
winds. 

(4) After the most massive stars evolved to the WR stage, charac- 
terised by dense and fast stellar winds, the level of cluster wind 
diffuse X-ray emission rises dramatically. Powered further by su- 
pernova explosions, the diffuse emission of a cluster wind is nearly 
constant during about 40 Myr. 

(5) After about 40 Myr of cluster evolution, the supernova rate is 
expected to drop steeply. At this stage, the cluster wind ceases. 



(6) The diffuse X-ray emission observed in the young Galactic clus- 
ter NGC 3603 cannot be explained by a c luster wind alone . Base d 
on a study of the Orion Nebular Cluster bv lFlaccomio et ai1l2003h . 
we model the evolution of X-ray emission from an ensemble of 
low-mass stars. We conclude that X-ray active low-mass stars are 
the dominant sources of X-rays in c lusters younger than » 2 Myr. 
Thus, we confirm the suggestion of iMoffat et alj l2002h that non- 
resolved low-mass stars are responsible for the rather high level of 
diffuse X-rays observed in NGC 3603. 

(7) We analyse the pointing Rosat and available XMM-Newton and 
Chandra observations of O and WN-type binaries and find that the 
correlation Lx = 10~ 'Lboi holds for these stars. 

(8) From the study of X-ray observations of WN8 type stars we 
conclude that stars of this spectral type show no evidence of be- 
ing X-ray sources. We speculate that the wind driving mechanism 
might differ between WN8 and other types of O and WN stars. 

(7) Assuming evolution with constant bolometric luminosity, the 
collective X-ray luminosity of massive stars is constant over the 
first ~ 2 Myr. After the most massive stars become X-ray dim WR 
stars, the collective X-ray luminosity of massive stars declines fast. 

(8) The ratio of diffuse luminosity of the cluster wind to the X-ray 
luminosity of massive stars is lower for less massive clusters. 

(9) If a young supernova remnant is present in a cluster, it will dom- 
inate the cluster X-ray emission, however only for a short time of a 
few hundred years. The X-ray emission from a supernova remnant 
in the Sedov-Taylor phase is of the same order of magnitude as the 
cluster wind emission, but is expected to show a softer spectrum. 

(10) Since the X-ray emission from a star cluster sensitively de- 
pends on its evolutionary stage, it may be used to constrain the 
cluster age and stellar population. 
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